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Abstract 

q/ ' We study the existence and the exponential ergodicity of a general interacting parti- 

P^ ■ cle system, whose components are driven by independent diffusion processes with values 

r~| ! in an open subset of IR,'^, d > 1. The interaction occurs when a particle hits the bound- 

"jrt I a-i'y: it jumps to a position chosen with respect to a probability measure depending on 

the position of the whole system. 

Then we study the behavior of such a system when the number of particles goes 
to infinity. This leads us to an approximation method for the Yaglom limit of multi- 
^ ' dimensional diffusion processes with unbounded drift defined on an unbounded open 

C^ , set. While most of known results on such limits are obtained by spectral theory argu- 

ments and are concerned with existence and uniqueness problems, our approximation 
method allows us to get numerical values of quasi-stationary distributions, which find 
l/"^ ■ applications to many disciplines. We end the paper with numerical illustrations of our 

^D I approximation method for stochastic processes related to biological population models. 

Key words : diffusion process, interacting particle system, empirical process, quasi-stationary 
distribution, Yaglom limit. 
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1 Introduction 

Let D C R'^ be an open set with a regular boundary (see Hypothesis [1]). The first part of this 
paper is devoted to the study of interacting particle systems (X^,...,X^), whose components 
X* evolve in D as diffusion processes and jump when they hit the boundary dD. More 
precisely, let A^ > 2 be the number of particles in our system. Let us consider A^ independent 
d- dimensional Brownian motions B^,...,B^ and a jump measure J'^^^ : d{D^) \-^ Aii{D^), 
where A^i(D^) denotes the set of probability measures on D^ . We build the interacting 
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particle system {X^,...,X^) with values in D^ as follows. At the beginning, the particles X* 
evolve as independent diffusion processes with values in D defined by 

dX« = dBl + qf\xi'^)dt, X« e D, (1) 

where q^ is locally Lipschitz on D, such that the diffusion process doesn't explode in finite 
time. When a particle hits the boundary, say at time Ti, it jumps to a position chosen 
with respect to J'^^\Xl_^_,...,X^_). Then the particles evolve independently with respect 
to ([T]) until one of them hits the boundary and so on. In the whole study, we require the 
jumping particle to be attracted away from the boundary by the other ones during the jump 
(in the sense of Hypothesis [2] on JT'^^-' in Section 12. 2p . We emphasize the fact that the 
diffusion processes which drive the particles between the jumps can depend on the particles 
and their coefficients aren't necessarily bounded (see Hypothesis [T]). This construction is a 
generalization of the Fleming- Viot type model introduced in |5| for Brownian particles and in 
J20|for diffusion particles. Diffusions with jumps from the boundary have also been studied 



in |3|, with a continuity condition on JT"^^-* that isn't required in our case, and in 19|, where 
fine properties of a Brownian motion with rebirth have been established. 

In a first step, we show that the interacting particle system is well defined, which means 
that accumulation of jumps doesn't occur before the interacting particles system goes to 
infinity. Under additional conditions on q^ and D, we prove that the interacting particle 
system doesn't reach infinity in finite time almost surely. In a second step, we give suitable 
conditions ensuring the system to be exponentially ergodic. The whole study is made possible 
thanks to a coupling between {X^,...,X^) and a system of A^ independent 1-dimensional 
reflected diffusion processes. The coupling is built in Section 12.31 

Assume that D is bounded. For all A^ > 2, let ^T" be a jump measure and (q^ )i<i<N 
a family of drifts. Assume that the conditions for existence and ergodicity of the interacting 
process are fulfilled for all A^ > 2. Let M^ be its stationary distribution. We denote by 
X^ the associated empirical stationary distribution, which is defined by X'^ = j^Y^i^i^xi, 
where {xi,...,Xn) G -D^ is distributed following M^. Under some bound assumptions on 
(% )i<i<N,2<N (see Hypothesis H]) , we prove in Section [23] that the family of random mea- 
sures Af^ is uniformly tight. 

In Section [3l we study a particular case: q^ = q doesn't depend on i,N and 

jr(^)(xi,...,x^) = ^±^J2^^.^ ^i e dD. (2) 

It means that at each jump time, the jumping particle is sent to the position of a particle 
chosen uniformly between the A^ — 1 remaining ones. In this situation, we identify the limit of 
the family of empirical stationary distributions {X^)n>2- This leads us to an approximation 
method of limiting conditional distributions of diffusion processes absorbed at the boundary 
of an open set of R*^, studied by Cattiaux and Meleard in [7| and defined as follows. Let 
Uoo C R"^ be an open set and P°° be the law of the diffusion process defined by the SDK 

dX^ = dBt + VV{X^)dt, X^ e f/oo (3) 



and absorbed at the boundary dUoo- Here i? is a ci-diniensional Brownian motion and V G 
C^iUoo,^)- We denote by rg the absorption time of the diffusion process ([3]). As proved in 
[T], the hmiting conditional distribution 



hmP~(XrG.|t<ra) (4) 



exists and doesn't depend on x G f/oo, under suitable conditions which allow the drift VV^ 
and the set Uoo to not fulfill the conditions of Section [2] (see Hypothesis |5] in Section [3]). 
This probability is called the Yaglom limit associated with P°°. It is a quasi-stationary 
distribution for the diffusion process ([3]), which means that P^ i^t" ^ ^,^1^ < '^d) = t^oo for 



alH > 0. We refer to p, |23|, |25| and references therein for existence or uniqueness results on 
quasi-stationary distributions in other settings. 

Yaglom limits are an important tool in the theory of Markov processes with absorbing 
states, which are commonly used in stochastic models of biological populations, epidemics, 
chemical reactions and market dynamics (see the bibliography |29|, Applications]). Indeed, 
while the long time behavior of a recurrent Markov process is well described by its stationary 
distribution, the stationary distribution of an absorbed Markov process is concentrated on 
the absorbing states, which is of poor interest. In contrast, the limiting distribution of the 
process conditioned to not being absorbed when it is observed can explain some complex 
behavior, as the mortality plateau at advanced ages (see jl| and (32|), which leads to new 
applications of Markov processes with absorbing states in biology (see |24|). As stressed in 
J28| . such distributions are in most cases not explicitly computable. In J7|, the existence of 
the Yaglom limit is proved by spectral theory arguments, which doesn't allow us to get its 
explicit value. The main motivation of Section [3] is to prove an approximation method of z/qo, 
even when the drift W and the domain Uoo don't fulfill the conditions of Section [2l 

The approximation method is based on a sequence of interacting particle systems defined 
with the jump measures ([2]), for all N >2. In the case of a Brownian motion absorbed at the 
boundary of a bounded open set {i.e. q = 0), Burdzy et al. conjectured in |j| that the unique 
limiting measure of the sequence {X'^)N^^i& the Yaglom limit i^oo- This has been confirmed 
in the Brownian motion case (see (SJ, (iSf and J26|) and proved in jl6| for some Markov 
processes defined on discrete spaces. New difficulties arise from our case. For instance, the 
interacting particle process introduced above isn't necessarily well defined, since it doesn't 
fulfill the conditions of Section [2J To avoid this difficulty, we introduce a cut-off of Uoo near 
its boundary. More precisely, let {Um)m>o be an increasing family of regular bounded subsets 
of f/oo, such that W is bounded on each Uoo and such that Uoo = Um>o ^o^- We define an 
interacting particle process (X'^'^,...,X'^'-'^) on each subset U^, by setting ql = W and 
D = Umin 1^. For all m > and iV > 2, (X'"'\...,X'"'^) is well defined and exponentially 
ergodic. Denoting by X"^'^ its empirical stationary distribution, we prove that 

lim lim A""'^ = z/oo. 

We conclude in Section 13.31 with some numerical illustrations of our method applied to 
the 1-dimensional Wright-Fisher diffusion conditioned to be absorbed at 0, to the Logistic 
Feller diffusion and to the 2-dimensional stochastic Lotka-Volterra diffusion. 



2 A general interacting particle process with jumps from 
the boundary 

2.1 Construction of the interacting process 

Let D be an open subset of R'^, d>\. Let A^ > 2 be fixed. For all i G {!,... ,A^}, we denote 
by P* the law of the diffusion process X^*\ which is defined on D by 

dXf = dBl - qf\xP)dt, X^'^ =x'eD (5) 

and is absorbed at the boundary dD. Here B^,...,B^ are A^ independent (i-dimensional 
Brownian motions and ql = {ql^ ,...,ql^ ) is locally Lipschitz. We assume that the process 
is absorbed in finite time almost surely and that it doesn't explode to infinity in finite time 
almost surely. 

The infinitesimal generator associated with the diffusion process (^ will be denoted by 
£(^\ with 



rW _ 1 V^_^ _ ( 






on its domain V (n). 
'-i 

For each i G {1,...,A^}, we set 

T^i = {{xi,...,xn) G d{D^), such that Xj G dD, and, Vj 7^ i, Xj G D}. 

We define a system of particles {X^,...,X^) with values in D'^ , which is cadlag and whose 
components jump from |J^ Dj. Between the jumps, each particle evolves independently of the 
other ones with respect to P*. 

Let J'^'^^ : Uj^o^i ~^ -^i(-D) be the jump measure, which associates a probability 
measure J'^'^\xi,...,Xn) on D to each point {xi,...,Xn) G IJj=i^i- Let {Xq,...,X^) G -D^ be 
the starting point of the interacting particle process {X^,...,X^), which is built as follows: 

• Each particle evolves following the SDE ([5]) independently of the other ones, until one 
particle, say X^^, hits the boundary at a time which is denoted by Ti. On the one hand, 
we have Ti > almost surely, because each particle starts in D. On the other hand, 
the particle which hits the boundary at time Ti is unique, because the particles evolves 
as independent Ito's diffusion processes in D. It follows that (X^^_,...,X^_) belongs to 

• The position of X*^ at time Ti is then chosen with respect to the probability measure 

j(^\xl,...,x^j. 

• At time Ti and after proceeding to the jump, all the particles are in D. Then the 
particles evolve with respect to ([5]) and independently of each other, until one of them, 
say X*2, hits the boundary, at a time which is denoted by T2. As above, we have ri < T2 
and(X4_,...,X,^_)Gl).,. 



• The position of X^^ at time T2 is then chosen with respect to the probabihty measure 

• Then the particles evolve with law P* and independently of each other, and so on. 

The law of the interacting particle process with initial distribution m G Aii{D^) will be 
denoted by P^, or by P^ iim = S^, with x G -D^. The associated expectation will be denoted 
by E^, or by E^ ii m = 6^. For all /3 > 0, we denote by Sfs = inf{t > 0, \\{Xi,...,Xn)\\ > 13} 
the first exit time from {x G D^ , \\x\\ < f3}. We set 6*00 = lim/3_5.oo ^/j. 

The sequence of successive jumping particles is denoted by {in)n>i, and 

< Ti <T2 < ... 

denotes the strictly increasing sequence of jumping times (which is well defined for all n > 
since the process is supposed to be absorbed in finite time almost surely). Thanks to the 
non-explosion assumption on each P*, we have r„ < S^o for all n > 1 almost surely. We 
set Too = lim„^oo Tn < 'S'oo- The process described above isn't necessarily well defined for all 
t G [0,5*00 [, and we need more assumptions on D and on the jump measure J'^'^' to conclude 
that Too = Soo almost surely. 

In the sequel, we denote by (po the Euclidean distance to the boundary dD: 

4>d{x) = inf \\y — x\\2, for all x & D. 

y£dD 

For all r > 0, we define the collection of open subsets Dr = {x E D, (J)d{x) > r}. For all 
/3 > 0, we set i?/? = {a; G -D, ||a;|| < /?}. 

Hypothesis 1. There exists a neighborhood U of dD such that 

1. the distance 0^ is of class C^ on U, 

2. for all /3 > 0, 

inf C- (J)d(x) > —00. 

x&UnBfi, i£{l,...,N} 

In particular. Hypothesis [1] implies 

\\VMx)\\2 = i,^xeU. (6) 

Remark 1. For example, the first part of Hypothesis [1] is fulfilled if D is an open set whose 
boundary is of class C^ (see jl2|. Theorem 4.3]). It is also satisfied by the rectangle with 
rounded corner defined in Section [3.3.31 

The following assumption ensures that the jumping particle is attracted away from the 
boundary by the other ones. 



Hypothesis 2. There exists a non- decreasing continuous function f^'^' : R+ — )■ IR+ vanishing 
at and strictly increasing in a neighborhood of such that, Vz G {1,...,A^}, 

^ inf :r(^)(xi,...,x^)({y G A 0z.(y) >min/(^)(0z.(x,))}) >pr\ 

Po > zs a positive constant. 

Informally, f^^\(t>D) is a kind of distance from the boundary and we assume that at each 
jump time r„, the probability of the event "the jump position X^^ is chosen farther from the 

boundary than at least one another particle" is bounded below by a positive constant p^ . 

Remark 2. Hypothesis [2] is very general and allows a lot of choices for J'^^\xi,...,Xn)- For 
instance, for all fi G A4i{D), one can find a compact set f^ C -D such that f^t{K) > 0. Then 
J^^\xu...,xn) = /u fulfills the assumption with pj,^^ = /i(K) and f^^\(t)D) =^d r\d{K,dD). 
Hypothesis [2] also includes the case studied by Grigorescu and Kang in |20|, where 

J^^\xi,...,xn) = ^Pij{xi)S^^, \/{xi,...,xn) G Vi. 

with Yljj=iPiji^i) — 1 ^^^ ^^Ue{i,...,N},jjti^XiedD Pij{xi) > 0. In that case, the particle on the 
boundary jumps to one of the other ones, with positive weights. It yields that Hypothesis [2] 
is fulfilled with Pq = 1 and f^^\(t)D) = 4>d- In Section [3l we will focus on the particular 
case 

J'^^\xi,...,Xn) = J^—r ^ 4^., V(Xi,...,XAr) G Pi. 

j=l,...,N, j^i 

That will lead us to an approximation method of the Yaglom limit (jlj). 

Finally, given a jump measure J'^^^ satisfying Hypothesis 2 (with Pq and f^^^), any 
0"*-^'' : IJi=o -^^ ~^ ■Mi{D) and a constant a^^^ > 0, the jump measure 

J-i^)(xi,...,x;v) = a(^):r(^)(xi,...,x^) + (1 - a(^))a(^)(xi,...,x^), V(a;i,...,a:^) G V^, 

fulfills the Hypothesis [2] with p'^^^J = a(^)pj^^ and f^^\(pD) = f^^\<pD)- 

Finally, we give a condition which ensures the exponential ergodicity of the process. In 
particular, this condition is satisfied if D is bounded and fulfills Hypothesis [TJ 

Hypothesis 3. There exists a > 0, tg > and a compact set Kq C D such that 

1. the distance (p^ is of class C^ on D\ D20 O'l^'d 

inf C- (pDix) > —00. 

x&D\D2c, ie{l,...,N} 



2. for all i G {!,..., A^}, we have 

N 

\ 

■ x£D^/2 






i=l 
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Theorem 2.1. Assume that Hypotheses IJ\ andlM are fulfilled. Then the process {X^,...,X^) 
is well defined, which means that Too = Soo almost surely. 

If Hypothesis [H and the first point of Hypothesis are fulfilled, then Too = Soo = +00 
almost surely. 

If Hypotheses\^ and\^ are fulfilled, then the process (X^,...,X^) is exponentially ergodic, 
which means that there exists a probability measure M^ on D^ such that, 

||PJ^((X^...,xf ) G .) - M^IItv < C(^n^) (P^^^^)* , Vx e D^, Vt e R+, 

where C^^\x) is finite, p^^'^ < 1 and \\-\\tv is the total variation norm. In particular, M^ 
is a stationary measure for the process (X^,...,X^). 

The main tool of the proof is a coupling between (X/,...,X/^)tg[o^5„] and a system of A^ 
independent one-dimensional diffusion processes {Y^' ,...,¥{' )t6[o,5o], for each (3 > 0. The 
system is built in order to satisfy 

< Yf'' < (pniXl) a.s. 

for all t G [0,Too A Sjs] and each i G {1,...,A^}. We build this coupling in Subsection 12.21 and 
we conclude the proof of Theorem 12.11 in Subsection 12.31 . 

In Subsection 12. 4[ we assume that D is bounded and that, for all N > 2, we're given JT''-^-* 
and a family of drifts (g| )i<i<Ar, such that Hypotheses [1], [2] and [3] are fulfilled. Moreover, 
we assume that a in Hypothesis [3] doesn't depend on A^. Under some suitable bounds on 
the family (g| )i<j<Ar, Ar>2; we prove that the family of empirical distributions {X^)j\f>2 
is uniformly tight. It means that, Ve > 0, there exists a compact set K (Z D such that 
E{X^{D \ K)) < e for all A^ > 2. In particular, this implies that {X^)n>2 is weakly 



compact, thanks to [22|. Let us recall that a sequence of random measures {'~fN)N on D 
converges weakly to a random measure 7 on D, if E{'-fi\f{f)) converges to E{'~f{f)) for all 
continuous bounded functions / : D — > R. This property will be crucial in Section [31 

2.2 Coupling's construction 

Proposition 2.2. Assume that HypothesisUl is fulfilled and fix 13 > 0. Then there exists 
a > 0, a N -dimensional Brownian motion {W^,...y/^) and positive constants Qi,...,Qn 
such that, for each i G {!,... JV}, the reflected diffusion process with values in [0,a] defined by 
the reflection equation (cf. ImJ) 

Yf'' = Yt + Wl- Qa + Lf - ir, Yi'' = min{a,MX^)) (7) 

satisfies 

0<Yf'' <(j)D{Xi)Aa a.s. (8) 

for all t G [0,roo A Si^[ (see FigureU\)- In dZD, -^*'° (resp. V'"-) denotes the local time ofY^'^ 
at {0} (resp. {a}). 




Figure 1: The particle X^ and its coupled reflected diffusion process Y^ 



Remark 3. If the flrst part of Hypothesis [3] is fulfllled, then the proof remains valid with 
(3 = oo and a = a (where a > is deflned in Hypothesis [3]) . This leads us to a coupling 
between X^ and Y°°'\ valid for all t E [0,Too A Soo[= [0,Too[. 



Proof of Proposition \2.2\ : The set Bp\U is a compact subset of D, then there exists a > 
such that Bi3\U G D2a- In particular, we have Bp \ D2a C f/, so that 0d is of class C^ in 

Bp \ D2a. 

Fix i G {!,... ,A^}. We deflne a sequence of stopping times {O^n such that XI G 5^ \ D2a 
for all t G [6'2„,6'2„+i[ and XI G Da for all t G [02n+i^(^2n+2[- More precisely, we set (see Figure 

m 

ei = inf {t G [0, + oo[, Xl G B„ \ DJ A r^ A ^^, 
e\ = inf {t G [to, + oo[, Xl G D^} A Too A Sp, 



and, for n > 1, 



^lxieBp\Da]^T^^S, 



^^„+i = inf {t G [4n, + oo[, Xi G D^} A Too A 5^. 



/3' 



The sequence {91^ is non-decreasing and goes to Too A Sp almost surely. 

Let 7* be a 1-dimensional Brownian motion independent of the process (X^,...,X^) and 
of the Brownian motion {B^,...,B'^). We set 




Figure 2: Definition of tlie sequence of stopping times (6'Jj),„ 



>o 



and, for all n > 0, 

^l = ^k .. + (^* - ^k J f°^ ^ ^ [^2n+l,^2n+2[, 

where j„i 'W(I)d{XI_) ■ dB\ has the law of a Brownian motion between times 9\^ and ^2n+i5 

2n 

thanks to (El). The process (Vr\...,Vr^) is yet defined for all t G [0,roo A Sf}[. We set 
^l = ^Las,- + ill - iLas,) for t G [Too A S,, + oo[ 

It is immediate that {W^,...,W^) is a A^-dimensional Brownian motion. 

Fix i e {!,..., A^}. Thanks to Hypothesis O there exists Q\ >0 such that 

inf CfUn{x)>-Qr. 

Let us prove that the reflected diffusion process V^'' defined by ([7]) fulfills inequality ([8]) for 
alHe [Q,T^ASp[. 

We set C = inf <0 < t < t^o A Sfs, Yf''^ > (j)j;){Xl) > and we work conditionally to (^ < 

Too A 5^. By right continuity of the two processes, 

< (J)d{XI) < if '^ < a a.s. 

One can find a stopping time (' g]C,Too A Sj^l, such that X* doesn't jump between ( and (' 
and such that F/'* > and X^ E Bp\ D2a for all t E [CC] almost surely. 
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Thanks to the regularity of </>/) on Bp \-D2a, we can apply Ito's formula to ((/)£)(Xj*))ig[(^^(^/], 
and we get, for all stopping time t G [C;C']) 

MXl) = Mxi) + J VMXl) -dBl + J cfUniXDds. 

But C and (' lie between an entry time of X* to Bf^ \ Da and the following entry time to D2a- 
It yields that there exists n > such that [C,C] C [6*2^)^271+1 [■ ^^ deduce that 

Mxi) - Yf'' = MXD - if '^ + 1* {c^Ud{xi) + Qf^)ds - Lf + Lif + Vf - q\ 



where Cf^^cpo^X'^ + Q^^ > 0, {LY)s>o is increasing and Vf = LJ'°, since V^^' doesn't hit 
between times ( and t. It follows that, for all t G [C)C']) 

Mxi)-Yf''^>Mxi)~Yl''^ 

>MXl-)-Y^^j^>Q. 

where the second inequality comes from the positivity of the jumps of (f)j:,{X^) and from 
the left continuity of F^'*, while the third inequality is due to the definition of C^. Then 
0£)(X*) — y''* stays non-negative between times C, and C,' ., what contradicts the definition of 
C,. Finally, C, = Too /\Sp almost surely, which means that the coupling inequality ([8]) remains 
true for all t G [0,roo ^Sj3[. D 

2.3 Proof of Theorem [2J] 

Proof that (X^,...,X^) is well defined under Hypotheses\^ and\^ Let X > 2 be the size of 
the interacting particle system and fix arbitrarily its starting point x G D^ . Thanks to the 
non explosiveness of each diffusion process P*, the interacting particle process can't escape 
to infinity in finite time after a finite number of jumps. It yields that Too < 5*00 almost surely. 

Fix /3 > such that x & Bp and define the event Cp = {too < Sp}. Assume that Cp 
occurs with positive probability. Conditionally to Cp, the total number of jumps is equal to 
+00 before the finite time Too- There is a finite number of particles, then at least one particle 
makes an infinite number of jumps before Too- We denote it by Zq (which is a random index). 

For each jumping time r^, we denote by cr^° the next jumping time of zq, with r„ < a^' < 
Too- Conditionally to Cp, we get a^° — r„ — > when ri — )■ 00. For all C^ function / with 
compact support in ]0,2a[, the process f [cf) £){X^^)) is a continuous diffusion process with 
bounded coefficients between r„ and (T^n~-> then 

sup |/(0b(X,^°))| = sup \f{MXT))-f{MX\))\ ^0, a.s. 

Since the process 0d(X*'') is continuous between r^ and cr^° — , we conclude that 0d(X^°) 
doesn't lie above the support of /, for n big enough almost surely. But the support of / can 

10 



be chosen arbitrarily close to 0, it yields that </)£)(X*°) goes to almost surely conditionally 

to C/3. 

Let us denote by (T^°)n the sequence of jumping times of the particle io. We denote by 
An the event 



where f'^^^ is the function of Hypothesis [2] . We have, for all 1 < A; < /, 






n=k 



where, by definition of the jump mechanism of the interacting particle system. 



<i-pr, 



by Hypothesis [2l By induction on /, we get 



I 



P\{\K\ <(l-p^)'"^Vl<A:</. 

\n=k / 

Since Pq > 0, it yields that 

^ f u n -4-) = »■ 

\fc>l n=fc / 

It means that, for infinitely many jumps r„ almost surely, one can find a particle j such that 
/ (0D(-^r„)) — ^b(^tII)- Because there is only a finite number of other particles, one can 
find a particle, say jo (which is a random variable), such that 

/(^)(0Z5(^^:)) < 0d(X;o), for infinitely many n > 1. 

In particular, lim^^oo (0d(^;^),/^^H0d(X^:;))) = (0,0) almost surely. But {f^'^Y^ is well 
defined and continuous near 0, then 

hm (0z5(X;«),0z5(X^«)) = (0,0) a.s. 



Using the coupling inequality of Proposition I2.2[ we deduce that 



C^c(lim(F/'^«,y/'^«) = (0,0)1. 
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Then, conditionally to Cp, y/^-^o g^j^^j y/3jo qj-q independent reflected diffusion processes with 
bounded drift, which hit at the same time. This occurs for two independent reflected 
Brownian motions with probability 0, and then for Y^''^° and Y^'^° too, by the Girsanov's 
Theorem. That implies Px{Cp) = 0. 

We have proved that Too > Sf^ almost surely for all /3 > 0, which leads to t^q > S^o almost 
surely. Finally, we get Too = >S'oo almost surely. 

If the first part of Hypothesis [3] is fulfilled, one can defined the coupled reflected diffusion 
Y°°'\ which fulfills inequality ([8]) with a = a and for all t G [0,roo A Soo[= [0,roo[. Then the 
same proof leads to 



{Too < +00} c ( hm (y~'^o,y-'^o) = (0,0)) 



Finally, we deduce that Too = 00 almost surely. D 

Remark 4. One could wonder if the previous coupling argument can be generalized, replac- 
ing (|5]) by uniformly elliptic diffusion processes. In fact, such arguments lead to the definition 
of y* as the reflected diffusion Y^ = J (f){Xl.)dW^ — Qit + L° — L", where is a regular func- 
tion. In our case of a drifted Brownian motion, is equal to 1 and Y^ is a reflected drifted 
Brownian motion independent of the others particles. But in the general case, the Y^ are 
general orthogonal semi-martingales. It yields that the generalization of the previous proof 
reduces to the following hard problem (see J3l|, Question 2, page 217] and references therein): 
"Which are the two-dimensional continuous semi- martingales for which the one point sets 
are polar ?". Since this question has no general answer, it seems that the previous proof 
doesn't generalize immediately to general uniformly elliptic diffusion processes. 

We emphasize the fact that the proof of the exponential ergodicity can be generalized 
(as soon as Tqo = S^o = +C)0 is proved), using the fact that (Y^^,...,Y/^)t>o is a time changed 



Brownian motion with drift and reflection (see 3l|, Theorem 1.9 (Knight)]). This time change 
argument has been developed in 20|, with a different coupling construction. This change of 



time can also be used in order to generalize Theorem 2.3 below, as soon as the exponential 
ergodicity is proved. 

Proof of the exponential ergodicity. It is sufficient to prove that there exists n > 1, e > and 
a non-trivial probability ^ on D'^ such that 

P,((X^^(^),...,X^(^)) eA)> e^{A), Vx e Ko, A e i3(D^), (9) 

with Kq = f Kq j , where t^ and Kq are defined in Hypothesis [3], and such that 

sup E^^k"') < 00, (10) 



where k is a positive constant and r' = minjra > 1, (X"^ (jv),...,X^(iv))ng]N £ Kq\ is the 
return time to Kq of the Markov chain (X^ ,j^.,...,X^,j^.)n)z]^. Indeed, Down, Meyn and 
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^,m)nm 



Tweedie proved in 13|, Theorem 2.1 p. 1673] that if the Markov chain {X^ (^),...,X^ 

is aperiodic (which is obvious in our case) and fulfills ([9]) and flTOj) . then it is geometrically 
ergodic. But, thanks to |l3|, Theorem 5.3 p. 1681], the geometric ergodicity of this Markov 
chain is a sufficient condition for {X^,...,X^) to be exponentially ergodic. 



We assume without loss of generality that Kq C Da/2 (where a is defined in Hypothesis 



[3]). Let us set 



'-0 



''0 

Thanks to Hypothesis [3l ?9 is a non-trivial probability measure. Moreover, iQ is clearly 
fulfilled with n = 1 and e = Hiiinf^eD,, P^(^?l) e i^J^^). 

Let us prove that 3^; > such that flTU|) holds. One can define the A^-dimensional 
diffusion {Y°°'^,...,Y°°'^) reflected on {0,a} and coupled with {X^,...,X^), so that inequality 
(j8]) is fulflUed for all t G [0, + oo[ and a = a. For all xq G -D^, we have by the Markov 
property 



:',...,x'',^,) G O > p.o(^>) G Da/2,yi) inf F,{xy ^ /rl^V 






TV 



> 



p,,(x>) G D«/2,vo n inf Px(^;(iv) G /rj- 



(JV)^ 



. 2:6-00/2 
1=1 



>p,„(x;,^) gd,/2,v^k 

where p-[ > is deflned in Hypothesis [31 It yields that 



(TV) 






p.„((xV)'-,^.V)) e O >ri^^p.o(0D(x;(,,) > «/2,vo 



i=l 

thanks to Proposition 12.21 A comparison argument shows that Pyoo,i(y^'* > a/2) > 

ig 



Po(F°°'^ > a/2). Then 



N 



i=l 



thanks to the strict positivity of the density of the law of F^'*, for all i G {1,...,A^}. Using 
the Markov property, we get, Vn > 1, 



*0 



P{r' > 2nC) > (1 - inf P.„((X\,^„...,X^,,,) G K^))P{r' > 2{n - 1)^^^; 
> (1 - inf P,,((Xi,,„...,X%) G 0)^ 
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where < mfxaeo^xoii^^ (n),---,^'^(n)) G Kq) < 1. It yields that there exists k > such 
that (^ is fulfilled. ° " D 



2.4 Uniform tightness of the empirical stationary distributions 

In this part, the open set D is supposed to be bounded. Assume that a jump measure J'^^^ 
and a family of drifts (q^ )i=i,...,Ar are given for each N > 2. 

Hypothesis 4. HypothesesUl andlE are fulfilled for each N > 2 and HypothesislE is fulfilled 
with the same a for each N > 2. Moreover, there exists r > 1 such that 



N 

sup — > r'^'*^* '' < +00, 



N>2 N ^ 



i=\ 



W_ .^.f _, _ ri^). 



where Ql ' = -mi^(,D\D^ Q (Pd{x). 

For all A^ > 2, we denote by m^ G Aii{D'^) the initial distribution and by fi'^{t,dx) 
the empirical distribution of the A^-particles process defined by the jump measure J'^^^ and 
the family (ql )ie{i,...,N}- Its stationary distribution is denoted by M^ and its empirical 
stationary distribution is denoted by X^: 



N 



i=\ 

where (x^,...,x^) is a random vector in D^ distributed following M^ . 

Theorem 2.3. Assume that Hypothesis [7] is fulfilled. For all sequence of measures m^ G 
Aiii^D^) and allt > 0, the family of random measures {^n^it^dx)) ^ is uniformly tight. In 
particular, the family of empirical stationary distributions (A*^) is uniformly tight. 

Proof. Let us consider the process (X^,...,X^) starting with a distribution m^ and its coupled 
process {Y°°'^^...X°°'^). For all t G [0,roo[, we denote by fi"^{t,dx) the empirical measure of 
{Yi°°' ,...,Y^°^' ). By the coupling inequality ([8]), we get 

/i^(t,D,^) < /i'^(t,[0,r]), Vr G [0,a]. 

Using the Markov property, we deduce that, for all s < t, 

Exi,...,xN {fi''{t - s,D^.)) < Ey^,.^^^^^y^,M (/i'^(t - s.D';)) a.s. 

Then, by a comparison argument, 

i?xi,...,xiv {li''{t - s,D^)) < Eo,...,o (/i'^(t - s,D':.)) a.s. 



1 ^ 



<j;^}_^Po{Y,.:<r)^.s. fill 
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Thanks to the Girsanov's Theorem, we have 



Qr-'l-s-iQl''')Ht-s)\ MqI''')'(^-^) 



\ t — s ' t — s t — s / 

where {w^,...,w'^) is a A^-diniensional Brownian motion. By the Cauchy Schwartz inequahty, 
we get 



PoO'S- < r) < Je„ {{s^_,.,,,,.^,,(P^^f) E„ ((.«r'»L.-<«'»'.'(.-.))^| 



- V^» (■'-L.+tiL-.-tjido.'-l) 



where the second inequahty occurs, since < 6 , ,ri,a _ri,o ([0,r]) < 1 almost surely and 

t — s' t — s t — s 

the process e^*^' "'t~2(Qi ) * jg the Doleans exponential of 2Q- wl, whose expectation is 1. 
Taking the expectation in (fTTl) , it yields that 

I 1 ^ 

E^. (/x^(t,D,^)) < JPo (5.j_^+L- -l;i([0'^]))]v E^'^''^"'^'^'"^^ WO<s<t. 



i=l 



Thanks to Hypothesis HJ there exists Sq G]0,t[ such that ;^ XI j=i ^ ^ '■'^' "* *■* ^°^ i^ uniformly 

bounded in A^ > 2. But Pq\5 t ,ji,a ,i,o ([0,rl) ) goes to when r — )■ 0, so that the 

family of random measures {fi^{t,dx))N>2 is uniformly tight. 

If we set m^ equal to the stationary distribution M^ , then we get by stationarity that 
X'^ is distributed as /i^(t,.), for all A^ > 2 and t > 0. Finally, the family of empirical 
stationary distributions {X^)n>2 is uniformly tight. D 

3 Yaglom limit's approximation 

We consider now the particular case J'^^\xi,...,Xj\f) = j^ J2k=i k^i^^k- ^^ ^^^^ jump time, 
the particle which hits the boundary jumps to the position of a particle chosen uniformly 
between the A^ — 1 remaining ones. We assume moreover that q\ = q doesn't depend on 
i,N. In this framework, we are able to identify the limiting distribution of the empirical 
stationary distribution sequence, when the number of particles tends to infinity. This leads 
us to an approximation method of the Yaglom limits (jl]), including cases where the drift of 
the diffusion process isn't bounded and where the boundary is neither regular nor bounded. 
Let Uoo be an open domain of R'^, with d >1. We denote by P°° the law of the diffusion 
process defined on Uoo by 

dX^ = dBt - VV{X^)dt, X^ = xe f/oo (12) 

and absorbed at the boundary dUoo- Here i? is a d-dimensional Brownian motion and V G 
C^(t/oo,lR). We assume that Hypothesis [S] below is fulfilled, so that the Yaglom limit 

z/oo = lim P^ (X~ G .|t < ra) , Vx G f/oo (13) 

i— >-+oo 
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exists and doesn't depend on x, as proved by Cattiaux and Meleard in [71, Theorem B.2]. We 
emphasize the fact that this hypothesis allows the drift W of the diffusion process flT2|) to 
be unbounded and the boundary dUoo to be neither of class C^ nor bounded. In particular, 
the results of the previous section aren't available in all generality for diffusion processes with 
law P°°. 

Hypothesis 5. We assume that 

1. F^{t9 < +oo) = 1, 

2. 3C > such that G{x) = \VV\'^{x) - AV{x) > -C > -oo, Vx G U^o, 

3. G{R) — !■ +00 as R ^ oo, where 

G{R) = inf {G{x); \x\ > R and x G f/oo} , 

4- There exists an increasing sequence {Um)m>o of bounded open subsets ofUoo, such that 
the boundary of Um is of class G^ for all m > 0, and such that IJm>o ^™ ~ ^'^■ 

5. There exists Rq > such that 

g-2y(a;)^^ < oo and 



L 



Uoan{d{x,dUo^)>Ro} 



/ (/ Pi°° {x, y) dy ] e~^^''Ux <oo. 



'Uoon{d{x,dUao)<Ro} 

Here p^°° is the transition density of the diffusion process flT2|) with respect to the 
Lebesgue measure. 

According to [Tj, the second point implies that the semi-group induced by P°° is ultra- 
contractive. The assumptions 1-4 imply that the generator associated with P°° has a purely 
discrete spectrum and that its minimal eigenvalue — Aqo is simple and negative. The last 
assumption ensures that the eigenfunction associated with — Aqo is integrable with respect to 
g-2y(a;)^^ Finally, Hypothesis [S] is sufficient for the existence of the Yaglom limit ( I13p . 

Remark 5. For example, it is proved in [7| that Hypothesis [S] is fulfilled by the Lotka-Volterra 
system studied numerically in Subsection 13.3.31 Up to a change of variable, this system is 
defined by the diffusion process with values in Uoo = R+, which satisfies 



2 8 8 2F/ 



t ~ "--^t 



'2 „ ^, /'\a2\3 



* * I 2 8 8 2Yy 



(14) 



and is absorbed at dUoo- Here B^,B^ are two independent one- dimensional Brownian motions 
and the parameters of the diffusion process fulfill condition (13Q|) . 
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In order to define the interacting particle process of tlie previous section, we work witli 
diffusion processes defined on Um, m > 0. More precisely, for all ?7i > 0, we denote by P"* 
the law of the diffusion process defined on Um by 

rfXf- = dBt - qm{Xf-)dt, X^- =xeUm (15) 

and absorbed at the boundary dUm- Here S is a d-dimensional Brownian motion and 
Qm • Um I— >■ R is a continuous function. We denote by £*" the infinitesimal generator of 
the diffusion process with law P™. For all m > 0, the diffusion process with law P"^ clearly 
fulfills the conditions of Section H For all N > 2, let (X'"'\...,X'"'^) be the interact- 
ing particle process defined by the law P™ between the jumps and by the jump measure 
i7*-™'^''(xi,...,XAr) = jv^ Xlfc=i fcT^i^^fc- -^y Theorem 12. 11 this process is well defined and ex- 



ponentially ergodic. 

For all m > and all iV > 2, we denote by ji^'^{t,dx) the empirical distribution of 
(Xr'\...,Xr'^), by M™'^ the stationary distribution of (X™'i,...,X'"'^) and by X"^'^ the 
associated empirical stationary distribution. 

We are now able to state the main result of this section. 

Theorem 3.1. Assume that Hypothesisl5\ is satisfied and that qm = VV^lfj^ for all m. > 0. 
Then 

lim lim X""'^ = z/oo, 

m— >-cxD N—^oo 

in the weak topology of random measures, which means that, for all bounded continuous 
function f : Uoo H- R+, 

lim lim EiX"^'''if)) = u^if). 

m— >oo N—>oo 

In Section 13. Ij we fix ?Ti > and we prove that the sequence {X'^'^)n^2 converges to a 
deterministic probability Um when A^ goes to infinity. In particular, we prove that Um is the 
Yagloni limit associated with P'", which exists by [7|. In Section [3^ we conclude the proof, 
proceeding by a compactness/uniqueness argument: we prove that {i'm)m>o is a uniformly 
tight family and we show that each limiting probability of the family {i'm)m,>Q is equal to the 
Yagloni limit z/qo- The last Section IX^ is devoted to numerical illustrations of Theorem 13.11 

3.1 Convergence of {X^'^)n>2: when m > is fixed 

Proposition 3.2. Letm > be fixed and let qm '■ Um i— )■ R &e a continuous function. Assume 
that fi'^''^{0,dx) converges in the weak topology of random measure to a random probability 
measure fim with values in Aii[Um), when N — )■ oo. Then, for all T > 0, fi'^'^{T,dx) 
converges in the weak topology of random measure to P™ (X^- G .iXx G Um) when N goes to 
infinity. 

Moreover, if there exists Vm ^ -^i(f^m) such that 

Vm = lim P;: (Xf™ G .|Xf™ G Um) , V/i G Mi{Um), (16) 

then the sequence of empirical stationary distributions (A'™'^)7v>2 converges to Vm. in the 
weak topology of random measures when N goes to infinity. 
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Remark 6. In Proposition 13.21 Um is the Yagloni limit and the unique quasi- stationary 
distribution associated with P™. For instance, each of the two following conditions is sufficient 
for the existence of such a measure: 



1. If Qm = 1[7^V\^, by [71]. This is the case of Theorem 13.11 



2. If q.m belongs to C^'°(f/„) with a > 0, by [17 
Proof of Proposition \3.S[ We set 

where Af = Yl^=i '^T„<t denotes the number of jumps before time t. Intuitively, we introduce 
a loss of 1/A^ of the total mass at each jump, in order to approximate the distribution of the 
diffusion process f lT5|) without conditioning. We will come back to the study of /i"^'^ and the 
conditioned diffusion process by normalizing z/"^'^. 

From the Ito's formula applied to the semi-martingale ^^'^{t^ip) = ;^ X]j=i ^(^t"'*); 
where tjj G C^([/m,Il), we get 



^ E /^'"'"'(^n-,^), (17) 



N _ 

0<T„<t 



where /A'^''^{t,ip) is the continuous martingale 

N d ,t 

and Ai^'^{t,ilj) is the pure jump martingale 

N 



^EE/ ^i^T-w 



^E E i^i^T 



^ -"'"'^«-,^) 



Applying the Ito's formula to the semi- martingale u"^'^{t,ip), we deduce from (fT7|) that 

0<T„<t 



Where we have 



A'^ 



\-ir) - \-i^ 

But 



and 






^ N ) \ N ) N -I \ N 

Then 

That imphes 

0<T„<t ^ ^ 

It yields that, for all smooth functions \l/(t,x) vanishing at the boundary of Um-, 

z/'"'^(t,^(t,.)) - ^^'"■^(0,^(0,.)) = / y'^'''{s, ^1^ + c^m{s,.))ds 

Jo 9s 

where N''^'^{t,'$) is the continuous martingale 



4=1 J=l " ^ ' J 



and A/'-''^(t,\l/) is the pure jump martingale 



^E E (^)"""'"(*«.-'fy-i^^'"'"«-.*K--)))- 



j=l 0<T^<i 
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Let T > be fixed. For all 6 > 0, define ^'^(t,x) = P^_tP^f{x), where / e C^{U^) and 
(P™)s>o is the semigroup associated with P™ : P"^f{x) = ExifiX^"")). Then \1/^ vanishes 
on the boundary, is smooth, and fulfills 

■^^\s,x) + ^A¥{s,x) + qUx)V¥{s,x) = 0, 
thanks to Kolmogorov's equation (see |14| . Proposition 1.5 p. 9]). It yields that 

Since {^^^) " < 1 a-s., we get 



E{Af'^'^{T,¥y) <^||V^ 



5||2 

oo 



< 



T 



N ^{T -t + 5)M 



2 

CX) 



(19) 



where Cm > is a positive constant. The last inequality comes from J30l . Theorem 4.5] on 
gradient estimates in regular domains of R*^. The jumps of the martingale A^-''^(t,\l/'^) are 
smaller than -|:||\l/'^||oo, then 



E 



^ 1 

.0<r„<T 


fN-V 




V N , 


- JV2 1 ^ 1 oo^ 




■))} 




Lir<TV A^ ) \ 





Then 



i?(Ar^'^(vl/,T)2)<^||vl/||L<^||/||L 



Taking t = T and 5 = ^, we get from (HE]), ([I9]) and ([20]) that 



(20) 



E 



z/-.^(t,P?/)-i.-.^(0,P- /) 



T+-, 



< 



'CmT + 4 



N 



Assume that / vanishes at dUm, so that / belongs to the domain of C^. Then \\PTf — f\\oo < 

N I 



1 i|£™J||^ and we have 



^Je (|z/™'^(r,/) - z/™.^(o,p^/)n < J 



CmT + A 



N 



oo + ^l|/:"/lloo^^^=^o. (21^ 
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By assumption, the family of random probabilities (i^™'^(0,.))7v>2 = (yU*"'^(0,.))Ar>2 converges 
to Hm- We deduce from (12T|) that 

E {ly^'^^iTJ)) ^Eif^UPTf)) , (22) 

for all / G C^iUm) vanishing at boundary. But the family (//'"'^(T,.)) is uniformly tight 
by Theorem 12.31 . It yields from (122|) that its unique limiting distribution is yUm(P^.). In 
particular, 

Af->oo 

But /i,^(P™l(/^) never vanishes almost surely, so that 

The family of random probabilities {X"^'^)n>o is uniformly tight, by Theorem 12.31 Let 
A'"* be one of its limiting probabilities. By definition, there exists a strictly increasing 
map (yj : IN h-)> IN", such that X"^'"^^^^ converges in distribution to X"^ when A^ — )■ oo. By 
stationarity, A:'™''^(^) has the same law as yu'"''^*^^)(T,.), which converges in distribution to 
P;p(X^'" e .|X^'" eUm), thanks to ([23]). But P^„(X^'" G .IX^™ e t/^) converges almost 
surely to Um when T — )■ oo, by (IT6|l . We deduce from this that A'™' has the same law as Um- 
As a consequence, the unique limiting probability of the uniformly tight family {X"^)fyf is z/^, 
which allows us to conclude the proof of Proposition 13.21 D 



3.2 Convergence of the family {i'm)m>o 

Proposition 3.3. Assume that Hypothesis is fulfilled and that qm = Wl-u^. Then the 
sequence {i'rn)m.>o converges weakly to the Yaglom limit u^o when m — )■ oo. 

Remark 7. Since qm = Wljj^, the operator £"* is symmetric with respect to the mea- 
sure e~'^^^^^dx, but this isn't directly used in the proof of Proposition 13.31 We mainly use 
inequalities from J7| that are implied by the ultra-contract ivity of P°° and the third point of 
Hypothesis [5l However, it seems hard to generalize this last hypothesis and its implications 
to diffusions with non-gradient drifts. 

Proof of Proposition WTSi For all ??i > and ttt, = oo, it has been proved in [7| that —C^* 
has a simple eigenvalue Am > with minimal real part, where £*"* is the adjoint operator 
of £™. The corresponding normalized eigenfunction rjm is strictly positive on Um, belongs to 
C2([/m,R) and fulfills 

^'^*Vm = ->^mVm and / r]m{xfda{x) = 1, (24) 

Ju„, 

where 
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The Yaglom limit z/^ is given by 

rijnlu da 
dUm = -J- "^ , Vm > or m = oo. 

In order to prove that {vrn)m>o converges to Vooi we show that {\m)m>o converges to Aoo- 
Then we prove that {f]mXumd<^)m>o ^^ uniformly tight. We conclude by proving that every 
limiting point rjda is a nonzero measure proportional to rjooda. 



For all 777, > or m = oo, the eigenvalue Am, of — £™-* is given by (see for instance 34, 
chapter XI, part 8]) 

Am = inf (C 



^""^m.^m- 



where C^ {Um) is the vector space of infinitely different iable functions with compact support 
in Um and {f,g)^m ~ lu /('")5'('")'^'^('")- For all G C^ (f^oo), the support of (p belongs to 
Um for m big enough, then C^ (f^oo) = Um>o ^o° (^m) since the reverse inclusion is clear. 
Moreover, if G C^ (Um), then £°°*0(x) = £"**0(x) for all x G Um- Finally, 

Aoo=mf inf (/:™*0,0),m 

= lim \ Am- 

m>0 

Let us show that the family (?7mlt/m'^cr)m>o is uniformly tight. Fix an arbitrary positive 
constant e > and let us prove that one can find a compact set K^ C Uo^ which fulfills 

f emtu^da < e, Vm > 0. (25) 

Let -Ro be the positive constant of the fifth part of Hypothesis [5l For all compact set K, we 
have 

/ Vm^u^da = / rjmda + / Vmda. 

JUoo\K J {d{x,dU„^)>Ro}nUm\K J {d{x,dUm)<Ro}nUm\K 

From the proof of P, Proposition B.6], we have on the one hand 



/ Vrnda<JI e-^n-)dx, 

J {d{x,dUm)>Ro}nUm\K y J {d{x,dUoo)>Ro}nU^\K 

which is smaller than e/2 for a good choice of K, say K'^, since the integral at the right-hand 
side is finite by Hypothesis [5l On the other hand 

/ rjmda <e'^/^e^'^K ( p^°°{x,y)dy] dx, (26) 

J {d{x,dUm)<Ro}nU,„\K J{d{x,dU^)<Ro}nUac>\K \Ju^ J 

where k = supm>o ||'7me~^||oo < oo thanks to (7|, and Am < Aoo for all m > 0. But the 
integral on the right-hand side is well defined by Hypothesis [Sj then one can find a compact 
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set K'l such that ([26]) is bounded by e/2. We set K^ = K'^UK'J so that ([25]) is fulfilled. Since 
inequality ([25]) occurs for all e > 0, the family {rimda)m>o is uniformly tight. Moreover, rjmda 
has a density with respect to the Lebesgue measure, which is bounded by /te"^, uniformly in 
?Ti > 0. Then it is uniformly bounded on every compact set, so that every limiting distribution 
is absolutely continuous with respect to the Lebesgue measure. 

Let Tjda be a limiting measure of {r]md(j)m>o- For ^H ^ C^{Uoo,R), the support of (f) 
belongs to Um for m big enough, then 






iv, >C°°0) = lim {r]m, C 

= hm (£-*r/^,0)^,„ 

m— >-oo ' 

= lim -A^(r/„,0) 

Thanks to the elliptic regularity Theorem, t] is of class C^ and fulfills C°°*r] = — Aqo''?- But 
the eigenvalue X^o is simple, then rj is proportional to rj^o- Let /3 > be the non- negative 
constant such that rj = Prjoo- In particular, there exists an increasing function : IN i— )■ IN 
such that r]^[rn)do' converges weakly to firjooda. 

Let us prove that /3 is positive. For all compact subset K C t/oo, we have 



(5 (r/oo,e^lx>,.^ = lim (r/^M,li^e^)^^ 



<f>{m) 



- }^^^oo^^'^^^^)^'^^'^'t>'^^))^,4>{^) 



> -(l-sup(r/„„lc;„\^r/^) ) 
where k = sup„>o ||^me~^||oo < oo. For all m > and all i? > 0, 



(27) 



where G and G are defined in Hypothesis [5l Let us prove that (jlmiGrim) ^ m ^^ uniformly 
bounded in m > 0. For all x G t/„7, (1211) leads to 



Then 



1 1 

-G{x)rim{x) = Xmr]m{x) + -e^(^)A(r/^e-^)(x) 



{Vm,G7]m)^^^ = Xm {Vm,Vm)^^^ + ^ / Vmix)e ^^""^ A^rj^e ^){x)dx 

|Vr/^(x)e-^(")|2da; 
< Ai, 
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where the second equahty is a consequence of the Green's formula (see [2, Corollary 3.2.4]). 

But G{R) goes to +00 when R — )■ 00, then there exists i?i > such that -^.^ {Vm,^\x\>RiGrim) < 

i. Since k = sup^>Q ||?7me~^||oo < 00, we deduce from fl28|) that 



4 ju 

But one can find a compact subset Ki C f/00 such that J^ l{\x\<:Rj^}\K^dx < ^, then we 
have from ( !27|) 

/3(^o,lx).>^. 
It yields that /9 > and Proposition 13.31 follows. D 

3.3 Numerical simulations 

3.3.1 The Wright-Fisher case 

The Wright-Fisher with values in ]0,1[ conditioned to be absorbed at is the diffusion process 
driven by the SDE 



dZt = v/^t(l - Zt)dBt - Ztdt, Zo = z G]0,1[, 

and absorbed when it hits (1 is never reached). Huillet proved in (2l| that the Yagloni limit 
of this process exists and has the density 1 — 1x with respect to the Lebesgue measure. In 
order to apply Theorem 13. ![ we define P°° as the law of X°° = arccos(l — 2Z). Then P°° is 
the law of the diffusion process with values in U^ =]0,7r[, driven by the SDE 

1 _ 2 cos X°° 

dX?^ = dBt TT^dt, X^ = X G]0,7r[, 

* 2sinXf° ' ° ^ ' ^' 

absorbed when it hits (tt is never reached). One can easily check that this diffusion process 
fulfills Hypothesis [5l We denote by i/qo its Yaglom limit. 

For all m > 1, we define Um =]~7^~ ~[- Let P"^ and Um be as in Section [3l We proceed to 
the numerical simulation of the A^-inter acting particle system {X"^'^,...,X"^'^) with m = 1000 
and A^ = 1000. This leads us to the computation of ElX"^'^), which is an approximation of 
z/qo- After the change of variable Z. = 2cos(X.), we see on Figure [2] that the simulation is 
very close to the expected result (2 — 2x)dx, which shows the efficiency of the method. 

3.3.2 The logistic case 

The logistic Feller diffusion with values in ]0, + oo[ is defined by the stochastic differential 
equation 

dZt = y^tdBt + {rZt - cZ^)dt, Zq = z > 0, (29) 

and absorbed when it hits 0. Here S is a 1-dimensional Brownian motion and r,c are two 



positive constants. In order to use Theorem 13. H we make the change of variable X. = 2yZ.. 
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0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 
Figure 3: E(X'^''^) in the Wright-Fisher case 



This leads us to the study of the diffusion process with values in U^ 
absorbed at and satisfies the SDE 

1 ^ "V^oo ^f vcx)\3'' 



10, + oof, which is 



dxr 



dBt 



2X^ 



rXr ^ c{Xl 



dt, X^ = xG]0, + oo[. 



2 4 

We denote by P°° its law. Cattiaux et al. proved in |6| that Hypothesis [5] is fulfilled in 
this case. Then the Yagloni limit z/qo associated with P°° exists and one can apply Theorem 
13.11 with Um =] — ,'m'[ for all m > 1. For all iV > 2, we denote by P™ the law of the 
diffusion process restricted to Um and by X"^'^ the empirical stationary distribution of the 
iV-inter acting particle process associated with P™. 

We've proceeded to the numerical simulation of the interacting particle process for a large 
number of particles and a large value of m. This allows us to compute E{X'"^'^), which gives 
us a numerical approximation of Vooi thanks to Theorem 13. II 

In the numerical simulations below, we set m = 10000 and A^ = 10000. We compute 
E(X"^'^) for different values of the parameters r and c in (129|) . The results are graphically 
represented in Figure HI As it could be wanted for, greater is c, closer is the support of the 
QSD to 0. We thus numerically describe the impact of the linear and quadratic terms on the 
Yaglom limit. 



3.3.3 Stochastic Lotka-Volterra Model 

We apply our results to the stochastic Lotka-Volterra system with values in D 
in [7|, which is defined by the following stochastic differential system 



Ri studied 



dzl = V^dBl + [nZ] - c,^{Zlf - c^^zlz^) dt, 
dZ^ = V^dB^ + {r,Z^ - c^iZlZ^ - C22iZ^y) dt, 
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12 3 4 5 6 7 

Figure 4: £'(A'"^'^) for the diffusion process ( l29l) . with different values of r and c 



where (i?^,i?^) is a bi-dimensional Brownian motion. We are interested in the process ab- 
sorbed at dD. 

More precisely, we study the process X°° = (Y'^X'^) = (2 a/ 2'^^ /71, 2 1^^2/72), with values 
in f/00 = R^, which satisfies the SDE flT^ and is absorbed at dUoo- We denote its law by 
P°°. The coefficients are supposed to satisfy 

cii,C2i > 0, C1272 = C2i7i < and C11C22 - C12C21 > 0. (30) 

In |7|, this case was called the weak cooperative case and the authors proved that it is a suffi- 
cient condition for Hypothesis [5]to be fulfilled. Then the Yagloni limit z/qo = limi_5.+oo ^T (^°° £ 
is well defined and we are allowed to apply Theorem 13. II For each m > 1, we define Um as it 
is described on Figure O With this definition, it is clear that all conditions of Theorems 12.11 
and 13.11 are fulfilled. 

We choose m = 10000 and we simulate the long time behavior of the interacting particle 
process with A^ = 10000 particles for different values of C12 and C21. The values of the other 
parameters are ri = 1 = r2 = 1, Cn = C22 = 1, 7i = 72 = 1- The results are illustrated on 
Figure O One can observe that a greater value of the cooperating coefficients — C12 = — C21 
leads to a Yaglom limit whose support is further from the boundary and covers a smaller area. 
In other words, the more the two populations cooperate, the bigger the surviving populations 
are. 



\t < ra) 
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ci2 = C21 =-0,08 



C12 = C21 =-0,56 



Cl2 = C21 =-0.80 



Cl2 = C21=-0,88 



C12 = C21 =-0,96 



C12 = C2i =-0,996 



1 3,0 



Figure 6: Empirical stationary distribution of the interacting particle process for different 
values of C12 = C21 
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